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Abstract 

We study soliton solutions of a modified non-linear Schrodinger (MNLS) 
equation. Using an Ansatz for the time and azimuthal angle dependence 
previously considered in the studies of the spinning Q-balls, we construct 
multi-node solutions of MNLS as well as spinning generalisations. 
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I. INTRODUCTION 



Solitons and instantons [1] are classical, non-singular, finite energy solutions of non- 
linear field equations. There exist topological solitons and non-topological solitons. The 
former are characterised by a conserved topological charge which results (in most cases) 
from a spontaneous symmetry breaking of the theory. Examples are monopoles, vortices 
and domain walls. Non-topological solitons [2] have a conserved Noether charge which 
results -in contrast to the case of topological solitons - from a symmetry of the Lagrangian. 
Examples of these are so-called Q-balls [3] which are solutions of a complex scalar field 
theory involving a 6 potential. 

The topic of classical field equations' solutions that are rotating has gained a lot of 
interest in recent years. It is well known that the inclusion of gravity can lead to a number 
of rotating solutions such as the famous Kerr-Newman family of rotating black holes as well 
as globally regular gravitating solutions namely so-called boson stars [4]. These are spinning 
solutions of the coupled Einstein and Klein-Gordon equations. We refer to them as spinning 
rather than rotating since they are not rotating in the sense of classical mechanics, but have 
a time and azimuthal angle dependence of the form exp(iut + iNO). 

However, classical solutions in flat space that are spinning have only been constructed 
very recently [9] . These are Q-balls which were rotated in the sense that they have the time 
and azimuthal angle dependence chosen for boson stars. Solutions with k = 0,1,2 nodes 
of the scalar field and N — 1, 2, 3 have been constructed nuemrically. These solutions are 
of particular interest because since, as it has been shown in [10], the topological solitons in 
SU(2) gauge theory cannot rotate. 

In this paper, we study spinning solitons of a modified non-linear Schrodinger equation 
(MNLS) which arises, in the ontinuum limit, in the studies of solitons on discrete quadratic 
[5,7,6] and hexagonal lattices [8], respectively. These equations were studied recently for 
d > 2 space dimensions in a recent paper [11]. In the present paper, we reconsider in detail 
the case d = 2 adopting the viewpoint of [9]. Namely, we have managed to construct excited 
solutions for both non-spinning and spinning solitons. Physical quantities characterising the 
excited solutions are compared to those characterising the fundamental solutions. 

Our paper is organised as follows: in Section II we define the model, in section III we 
study its non-spinning solutions, while in section IV we concentrate our attention on the 
spinning generalisations. We present our conclusions in section V. 

II. MODIFIED NON-LINEAR SCHRODINGER EQUATION 

The modified nonlinear Schrodinger (MNLS) equation in 2 dimensions was presented in 
[5-8] and reads: 

+ A^ + ag^{\^\ 2 + b/\\^\ 2 ) =0 . (1) 

The constants a and b depend on the underlying geometry of the lattice, where a = 2, 
b = 1/12 for the quadratic lattice [5-7] and a = 4, b = 1/8 for the hexagonal lattice [8]. In 
the following, we will choose a, b according to the hexagonal lattice. 
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g is the coupling constant of the system and A denotes the Laplacian in 2 dimensions. 
The equation (1) can be derived from the functional defined by the following Lagrangian 
density: 

C = i (^V - ^*) + iWf - 2<#| 4 + ^(VH 2 ) 2 • (2) 

(In fact — £ is the Lagrangian density but we use the convention above throughout the 
paper.) The system has many conserved quantities [11], namely the space integral of the 
static part of the Lagrangian (2) which can be interpreted as the "static energy" : 

E = j ' (W| 2 - 2g\^ + \g{V\^\ 2 ) 2 ^ dxdy (3) 

as well as the norm n of the solution: 

n 2 = J ipip*dxdy . (4) 

We notice that by rescaling tp — > ip/^/g, we can scale out g from (1). (Of course, E and n 
are also rescaled then.) 



III. NON-SPINNING SOLITONS 
A. Ansatz 

ip is a complex scalar field, we choose first to consider it to be given by the (radially 
symmetric) Ansatz: 

^(t,r)=e^0(r) (5) 
where r = \/x 2 + y 2 . First, we remark that the quantity 



% J V dt dt , 
is conserved for the Ansatz (5), since 

Q = 2un 2 . (7) 

So, the action of the solutions constructed within the above Ansatz appears as a sum of two 
conserved quantities. 

Inserting the Ansatz (5) into (1), we obtain the following ordinary differential equation: 



„ | (j)' g<t><$>' 2 4g0 3 — uj(f) 
r 1 + g<p 2 1 + g<f) 2 



<f> + ~ + 1 I —2 + i I 172 = ' ( 8 ) 



where the prime denotes the derivative with respect to r. This equation can also be derived 
from the following effective one-dimensional Lagrangian density: 
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£ ns = (l + S0 2 )0 ,2 -2# 4 + cc>0 2 . (9) 

The form of the potential V{4>) = 2g(f) 4 — uocf 2 is given for g = 10 and u — 1, 0, —1 in 
Fig. 1. If we choose either g < or uo < 0, this potential has no local maxima. In a 
classical mechanics analog, we can think of this as a particle with position dependent mass 
in a x A potential. In order to have localised, finite energy solutions the function 0(r) has to 
approach a zero of the potential for r — > oo. This means here that <p(r — > oo) — > 0. If we 
come back to the analog with the particle, this means that (with r being now the "time") 
the particle starts of with zero velocity 4>'(0) = (this in fact is one of the used boundary 
conditions). It should end up at the local maximum of the potential (see Fig. 1). Employing 
the same argument as in Ref. [9]. the form of the potential suggests further that solutions 
with k = 1,2, ... nodes in the function <f>(r) should also exist for u > 0. 

It is convenient to also define the total energy of the solution according to: 

Etot = [ C ns rdrdd . (10) 



Our numerical results show that E tot is always positive (see the Numerical results section 
below) and E tot (k = 0) < E tot (k > 0). 

Before we discuss the numerical results we present the analytical behaviour of the solu- 
tions around the origin and at infinity. Close to the origin r = 0, the function 0(r) behaves 
like: 

0(r << 1) = C + 0(r 2 ) , C constant . (11) 

The asymptotic behaviour of the function <f>(r) can easily be determined from (8). The 
linearized equation is indeed a Bessel equation and leads to the following asymptotic be- 
haviour: 

0(r » 1) = J (x) ~ C ° S{X Z_ n/4) , if a; < (12) 

v 



<f>(r » 1) = K (x) ~ FV ^ ' , ifu>0 (13) 

where x = V^M 7 "- The function 0(r) therefore oscillates around = for to < and decays 
exponentially for oj > 0. Both type of behaviours are confirmed by our numerical analysis. 

Note that if the factor (1 + g<j) 2 ) is replaced by 1, the Lagrangian (9) reduces to that of a 
scalar field in a 4 potential. This is known to have no stable solutions, thus a 6 potential 
was introduced [3] in this case. In our model no such 6 term is possible, but we have a 
non-constant prefactor of the kinetic term, which leads to the existence of stable solutions. 
In the non-linear <7-model also a term in front of the kinetic term appears which results from 
the fact that the scalar field is constrainted. 

To construct soliton solutions of our model, we introduce the following boundary condi- 
tions for non-spinning solutions [8]: 

d r( j)(r = 0) = , 0(oo) = . (14) 

With these boundary conditions, the function 0(r) can cross the r-axis k times, i.e. multi- 
node solutions are possible. 
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B. Numerical results 

Without loosing generality, we can set g — 1 in the following. 

In Fig. 2, we present the two conserved quantities E, n and the value of at the origin, 
0(0), as function of uj for the non-spinning, zero- node soliton solution. We find that the 
norm for uj « 1 tends to a finite value as uj — > 0. It is interesting to connect this result 
with the one of [8]. With respect to this paper, we have rescaled ip — > ift/y/g, so that the 
norms are related by n = y/gn, where n denotes the norm of the unsealed solution. In [8], 
the convenion n — 1 was adopted, thus n = yfg and in the limit uj — > 0, the critical values 
of the norm is equal to \J 'g cr (N = 0, k = 0) ~ V2.93 ~ 1.71, which is exactly what we find. 

In [9], solutions were shown to exist only for a specific range of uj. Here, we find no upper 
or lower bound on uj as long as uj > 0. For uj < 0, the equation becomes of a Bessel type 
and the solutions become oscillatory. The difference between [9] and our results is due to 
the fact that in [9] a 6 -potential was introduced, while we have a non-constant prefactor in 
front of the kinetic term. The 6 -potential leads to the restriction in uj, while in our model 
the 4 -potential is not subject to such restrictions. 

The first excited solution has a maximum at r = 0. It then crosses zero at r = r and 
attains a minimum at , say r = r m . The numerical results show that (f>(r m ) ~ — 0(O)/2 and 
r m ~ 2r . For uj » 1 both r and r m decrease, so that the soliton becomes increasingly 
concentrated around its center. If uj « 1 these values get larger and the soliton is more 
delocalised. Since the norm of the first excited solution also tends to a finite value in the limit 
uj — > 0, we can conclude (following the above arguments for the fundamental solution and 
adopting the conventions of [8]) that the first excited solution exists only for a sufficiently 
high value of g, namely for y/g > \J g cr (N = 0, k = 1) ~ 4.38. This has to be compared to 
xj 'g cr (N = 0, k = 0) 1.71 of the fundamental solution. 

Finally, we present the total energy E tot of the solutions in Fig. 5. Clearly, the energy 
stays positive for all uj > and the zero-node (k = 0) solution has for all values of uj smaller 
energy than that of the first excited (k = 1) solution. 



IV. SPINNING SOLITONS 

To construct spinning generalisations, we take the following Ansatz [9] : 

rj, = (j)(r)e iwt+im , N integer . (15) 
The equation then takes the form: 

,// V aW 2 4# 3 -a>0 iV 2 

9 r 1 + gef) 2 1+gcj) 2 r 2 (l + gef) 2 ) ' 1 ' 

As before, this equation can equally be derived from the following static Lagrangian: 

£ , = i^V_ # < + ^ + ^. (17) 

where now the potential becomes explicitely r-dependent due to the additional term 
— iV 2 2 /(2r 2 ). Note that the equality (7) still holds for this Ansatz. The total energy 
E tot of the solution again is the space-integral of the static Lagrangian C s . 
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For r << 1, the function 0(r) behaves as: 

0(r « 1) = r N (C N + 0{r 2 )) , C N constant . (18) 

Moreover, the asymptotic behaviour of <f>(r) is now determined by the function Jn{x) and 
Kn(x) respectively for uj < and uj > 0. 

As a consequence (1) has to be solved subject to the following boundary conditions: 

0(r = 0) = , 0(oo) = . (19) 

Again, the choice of these conditions allows for multinode solutions. 

A. Numerical results 

Like for the non-spinning solutions, we choose g = 1. 

The conserved quantities E, n and the value of the derivative of 0(r) at the origin, 0'(O), 
for the N — 1 solutions with no nodes (k = 0) and one node (k — 1), respectively, are shown 
in Fig. 3. Again, the norm of the solutions tends to a finite value for uj — > 0. Thus, again 
employing the language of [8], we conclude that the normalised spinning solutions exist only 
for coupling constants g larger than a critical value g cr {N, k). Like in the non-spinning case, 
we find that the critical value of the k = solution is smaller than that of the k = 1 solution. 
In summary, we find for the critical values: 

g cr (0, 0) < Ml, 0) < g cr (0, 1) < 9cr(l, 1) • (20) 

The profiles of the solutions with N = 1 and k — 0, 1 are shown in Fig. 4 together with 
the two non-spinning solutions. The functions 0(r) have a maximum at r = r M and we find 
that the value of this maximum increases slowly as a function of uj. Again the solution and 
its first excitation gets narrower (resp. more spread out) for uj » 1 (resp. uj << 1). 

In Fig. 5, we present the total energy E tot of the solutions as a function of uj. Again, like 
in the = case, the energy is always positive and the energy of the fundamental solution 
with no nodes is smaller than that of the first excited solution. Comparing the energy of the 
solution for fixed k and different N, we find that the non-spinning solution has lower energy 
than the spinning solution, which is an expected result since the spinning adds energy to 
the solution. 

V. CONCLUSIONS 

The study of rotating solitons in flat space is still a little studied problem in classical 
field theory. Recently, non-topological spinning Q-balls have been constructed [9] using an 
Ansatz used to construct spinning boson stars [4]. 

In this paper, we have considered an equation which is related to the continuum limit 
of a system of equations describing the interaction of a complex Schrodinger-like field with 
a quadratic [5-7], resp. hexagonal lattice [8]. Although the system is originally not rela- 
tivistically invariant, the radially symmetric equations resemble classical static equations of 
a non-linear sigma model type supplemented by a "symmetry breaking" quartic potential. 
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The pattern of solutions of this non-linear equation seems extremely rich, containing 
several families of solutions (for this see also [11]). The solutions are mainly characterized 
by the parameter uo defining the harmonic time dependence. Remarkably, the function 
determining the soliton tends to zero in the limit oo — > although the "norm" of the solution 
stays finite in the same limit. This phenomenon provides a natural explanation of the fact 
that the solitons studied in [8] exist only for large enough values of the coupling constant. 

Among possible generalisation of our result we mention the construction of non-spinning 
and spinning solutions (and their excited solutions) of a three dimensional version of the 
equation and/or the coupling to a electromagnetic field, thus leading to charged solutions. 

Acknowledgments Y. B. gratefully acknowledges the Belgian F.N.R.S. for financial sup- 
port. B. H. was supported by an EPSRC grant. 
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FIG. 1. The potential V((f>) is shown for 5 = 10 and u = 1. 
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FIG. 2. The conserved quantities £7, n and the value (f)(0) are shown as functions of uj for the 

non-spinning solutions with k = 0, k = 1, respectively. 
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FIG. 4. The profiles <p(r) of the soliton solutions are shown for different choices of the node 
number k and the parameter N. N = are non-spinning solitons. g = 1. 
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